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We report specific heat and thermal conductivity of gadolinium- and yttrium-doped amorphous 
silicon thin films measured using silicon-nitride membrane-based microcalorimeters. Addition of 
gadolinium or yttrium to the amorphous silicon network reduces the thermal conductivity over a 
wide temperature range while significantly increasing the specific heat. This result indicates that a 
large number of non-propagating states are added to the vibrational spectrum that are most likely 
caused either by localized vibration of the dopant atom in a Si cage, or softening of the material 
forming the cage structures. High-resolution cross-sectional electron micrographs reveal columnar 
features in Gd-doped material which do not appear in pure amorphous silicon. Scattering from both 
the nanoscaled columns and the filled-cage structures play a role in the reduced thermal conductivity 
in the rare-earth doped amorphous semiconductor. The overall result is an amorphous solid with a 
large bump in C/T 3 and no plateau in thermal conductivity. 



PACS numbers: 61.43.Dq, 65.60.+a, 63.50.+X, 66.70.+f 
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I. INTRODUCTION 

A wide range of amorphous materials have similar fea- 
tures in specific heat, thermal conductivity and various 
spectroscopies that suggest a common physical origin. 
The most notable of these features are a linear term in 
the specific heat, C, usually observed below 2 K, thermal 
conductivity k ex T 1 - 8 in the same temperature range, 
a broad peak in C/T 3 between 10 and 50 K which is 
larger than in the corresponding crystalline phase, and 
a plateau in k over the same temperature rangei The 
peak in C/T 3 is often correlated with excess vibrational 
density of states seen in neutron or Raman scattering, re- 
ferred to either as a "boson peak"— or "excess modes. "^4 
Despite several decades of study, no microscopic the- 
ory exists which offers a complete explanation of these 
phenomena. The standard tunneling model is a widely 
accepted description of the behavior of C and k below 
2 K. In this model the linear term in C is attributed 
to a constant density of Schottky-like two-level systems 
(TLS), with the behavior of k explained by scattering 
from these states^ 6 - This model offers little explanation 
of the physical nature or origin of the very closely spaced 
energy levels. There have been several efforts to develop 
physical models which also explain the C/T 3 peak and 
k plateau. Examples are the soft-potential modelfi^ the 
fracton modeljSjifl and a model relating dynamics at the 
glass transition to the low temperature phenomena^i^ 

Amorphous silicon (a-Si) forms a tetrahcdrally bonded 
continuous random network and can be made only 
in thin-film form by evaporation, sputtering, and var- 
ious chemical vapor deposition techniques. Under- 



coordinated Si atoms form dangling bond defects, which 
can be reduced by introducing hydrogen either during 
deposition or after film growth. The hydrogenated mate- 
rial, a-Si:H, shows improved electronic properties and has 
several important industrial applications for large area 
microelectronic devices. a-Si is often studied theoreti- 
cally, in part because of its relatively simple structure, 
but also due to several predicted and observed devia- 
tions from typical amorphous solids. Both theory and 
experimental evidence suggest a strikingly low density of 
two-level systems, a relatively large change in the De- 
bye temperature between crystalline and amorphous sil- 
icon, and a relatively small peak in C/T 3 in a-Si, de- 
viating less from the Debye model than the crystalline 
phase. 4 - 5-13 ' 14 ' 15 Properties ranging from electrical con- 
ductivity to density of TLS also show a strong depen- 
dence on the a-Si him growth method. 16 Another exam- 
ple is the thermal conductivity plateau, though due to the 
difficulty of thermal measurements on thin- film samples, 
measurements of k in the expected temperature range of 
the plateau exist for only three samples: an w 50 /jm 
thick sputtered sample reported by Pompe et alii that 
shows a fairly well defined plateau, and much thinner 130 
and 277 nm thick films recently reported by our group 
grown by e-beam evaporation that show no evidence of 
the plateau^ Both results can be reasonably explained 
by the theory presented by Feldman 4 with different treat- 
ment of low-Q scattering. 

Amorphous rare-earth-silicon alloys (a-RE^Sii-j,, with 
RE=Gd, Y, Tb, etc.) prepared by e-beam co-deposition, 
have shown many phenomena related to the interac- 
tion between local RE magnetic moments and conduc- 
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tion electrons including a very large negative magne- 
toresistance at low temperatures! 1 ^* A negative Hall 
coefBcienljSii and thermopower 21 indicate that introduc- 
ing the large, heavy RE adds carriers to the electron 
band, with the Gd ions contributing localized S = 7/2 
magnetic moments, while Y provides a non-magnetic 
counterpart with nearly identical ionic radius and va- 
lency. Recent computational and experimental results 
shed some light on the structure of a-Gd x Sii_ x and its 
non-magnetic analog, a-Y^Sii-^. Local density func- 
tional theory simulations of a-Y^Sii-^ suggest that Y 3+ 
ions are surrounded by low-coordinated Si, leaving them 
in a "cage" of dangling bonds' X-ray absorption fine 
structure (XAFS) experiments on a-Gd x Sii_ x using both 
the Si and Gd absorption edges give a similar picture, 
with Si atoms as nearest-neighbors to the Gd ions and 
no clustering!* 3 - The XAFS studies also indicate that nei- 
ther the Si coordination number nor the Gd-Si distance 
changes with Gd composition, suggesting a Si cage but no 
dangling bonds. Electron spin resonance (ESR) measure- 
ments indicate that the addition of even small amounts 
of Gd to the Si matrix eliminates the dangling bonds. 24 
Cages in a-Gd x Sii_ x are further supported by calcula- 
tions using the full potential linearized augmented plane 
wave (FLAPW) method^ 

This filled-cage structure is reminiscent of the filled 
skutterudite antimonides, such as CeFe4Sbi2** or 
REi_ 2/ Fe 4 - :c Co a; Sbi2 (RE=La, Ce)*£, or Sr or Eu doped 
Ge clathratesjS* where the heavy La, Ce, Sr, or Eu 
dopant atoms fill cages in the host crystal structure. 
Slack originally proposed that "rattling" movement of 
the heavy dopant atom in the anharmonic potential of the 
cage would strongly scatter phonons but not charge car- 
riers, resulting in a material with a high thermoelectric 
figure of merit i2* Although there is some debate about 
the detailed mechanism of the phonon scattering) 3 *!*^* it 
is clear that the filled-cage structures in these materials 
dramatically alter the vibrational spectrum, reducing the 
thermal conductivity over a wide temperature range (at 
least 2 — 300 K), and increasing the specific heati 30 i 32 i 33 

In this paper we present specific heat and thermal 
conductivity measurements of gadolinium- and yttrium- 
doped amorphous silicon thin films from 3—100 K, as well 
as high-resolution cross-section transmission-electron mi- 
croscope (XTEM) observations of the films. We compare 
these measurements to data on a-Si films grown by the 
same technique and literature values for crystalline sili- 
con to probe the nature of vibrational states and scatter- 
ing in rare-earth doped amorphous silicon. 



II. EXPERIMENT 

Thin film a-Y x Sii— x and a-Gd^Sii-^ samples were 
e-beam co-evaporated from separate Y, Gd, and Si 
crucibles onto amorphous Si-N membrane-based mi- 
crocalorimeters and amorphous Si-N coated Si sub- 
strates. The microcalorimeters and substrates were held 



FIG. 1: Low Magnification XTEM images comparing a) a-Si 
and b) a-Gdi8Sis2 thin films. a-Si shows no evidence of crys- 
tallinity. a-Gdi8Sis2 also shows no evidence of crystallinity, 
but does reveal columnar structures with several nanometer 
length scale. 

near room temperature throughout the deposition, pro- 
moting the growth of amorphous films. Typical depo- 
sition pressures were < 2 x 10~ 8 Torr. The films on 
substrates were used to measure the film thickness via 
profilometry, the composition by Rutherford backscat- 
tering (RBS), the dangling bond density via ESR, the 
sound velocity^ and for XTEM imaging. Both XAFS 23 
and RBS measurements indicate atomic densities con- 
sistent with pore-free films. A picosecond ultrasonic 
measurements** of the longitudinal sound velocity in a a- 
GdzSii-z film gave v L = (5.39 ±0.44) x 10 5 cm/s, lower 
than the value measured in our a-Si, vl = (7.51 ±0.30) x 
10 5 cm/s. Detailed description of the microcalorimetry 
techniques for measuring C and k, including the deter- 
mination and subtraction of background contributions, 
appear elsewhere ,2****?*^ 

III. RESULTS AND DISCUSSION 

Figures and [21 compare XTEM images for a-Si and 
a-GdisSis2 films. The low magnification images (Fig. 
^) show a featureless a-Si film, whereas the a-GdisSis2 
film grown under identical conditions shows a vertically 
streaked appearance indicative of a columnar structure. 
The higher magnification images in Fig. [21 also show fea- 
tureless amorphous silicon, and 3 — 4 nm columnar fea- 
tures in amorphous gadolinium-silicon. Note that all im- 
ages suggest dense, pore-free, clearly amorphous films at 
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FIG. 2: High resolution XTEM images of the same films clar- 
ify the results of Fig. The approximate spacing of the 
columnar features in a-Gdi8Sis2 is 3 — 4 nm. 



the atomic level. Columnar microstructure, such as that 
seen in a-Gdi8Sis2, is a common outcome of the vapor 
deposition process for evaporated films when atomic mo- 
bility at the growth surface is low. 

Figure |2l compares k of a-Y^Sii-^ and a-Gd x Sh-x 
films to Pompe et al.'s sputtered 50 /im thick a-Si 
filmii and e-beam evaporated thin-film a-Si. The top 
axis indicates the estimated wavelength of the vibra- 
tions which carry the most heat in the dominant phonon 
approximation,-- calculated for a-Y^Sigg- The e-beam 
a-Si shows a lower k than the extremely thick sputtered 
film and no plateau. Addition of Y or Gd to the mate- 
rial further reduces k over the entire temperature range 
measured, and the three alloy films all have the same k 
within error bars at all T. This reduction occurs despite 
the addition of electrons to the material, which increases 
the electrical conductivity dramatically. However, this 
is not surprising, since a Wiedemann-Franz law estima- 
tion of the electronic contribution suggests that k is to- 
tally dominated by vibrational excitations^ In addition 
to the lack of composition dependence, k of a-Gd x S\i- x 
(which shows a huge negative magnetoresistance at low 
T) showed no measurable change in applied magnetic 
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FIG. 3: k vs. T plot comparing a-Yi2Siss, a-GdioSigo, and 
a-Gdi4Si86 to thin e-beam evaporated a-Si films^" and pre- 
viously published values for a-Si in the plateau regime^- 
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FIG. 4: C/T 3 vs. T plot showing a-Y^Sh-*, a-Sir^ and 
rr-Si4- The dashed line indicates the measured Debye contri- 
bution in a-Si, and the solid line through the a-Y x Sh-x data 
is a fit described below. 



fields up to 8 Tesla. 

Figure 0] compares C/T 3 vs. T (in J/g K 4 ) for a- 
YgSitji to thin-film a-Sii£ and bulk crystalline silicon. 
The dashed line is the Debye specific heat function, Cd, 
for 9 D = 487 K. C/T 3 for a-YgSigi shows a much larger 
bump than that seen in either a-Si or crystalline Si. a- 
Y^Sii-a; for a broad range of x have very similar specific 
heat above 60 K, and all show a large maxima in C/T 3 
equal to or greater than that shown here. Samples with 
larger x have electronic terms, 7T (which appear as 7/T 2 
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FIG. 5: k vs. T for dense, pore-free a-Si and a-Y^Sii-z thin 
films compared to data for Vycor, a silica glass with 29% 
porosity (dashed line) |— and a borosilicate glass sample with 
artificially structured arrays of pores^- k of a-SiC>2 is not 
shown but is only slightly higher than a-Si values at 10 K. 1 
At this temperature k for these two porous materials is 2 — 4 x 
lower than for a-Si02, while adding Y to a-Si reduces k by a 
factor of ~ 3. 



on a C/T 3 plot), and an additional contribution to C oc- 
curs in samples near the metal-insulator transition^ We 
have previously reported the specific heat of a-Gd^Sii-^, 
which is similar to values for a-V^Sii-^ above 60 K, but is 
dominated at lower temperatures by large contributions 
from magnetic degrees of freedomi22iii42i 

Figure compares k vs T below 100 K for a-Si and 
a-Y x Sii_ x to two porous glasses, a porous silica (Vycor) 
sample with 29% porosity^ (dashed line), and a glassy 
borosilicate material formed by fusing an array of capil- 
laries (V)ii The data for both porous materials has been 
corrected for the missing volume of the pores. Neither of 
these materials have a k plateau in the typical tempera- 
ture range, though Vycor shows an apparent plateau at 
much lower temperature. The reduction of k in these 
porous materials is due to enhanced damping or scatter- 
ing of long-wavelength modes caused by the poresM^ 
The addition of Gd- and Y- to a-Si reduces k by a quan- 
titatively similar factor as does the addition of pores to 
a-SiC>2, but must rely on a different mechanism for scat- 
tering of heat carriers. 

The structural evidence from XTEM (shown in Figs. 
Q] and 12), XAFS measurements, and simulations suggests 
two likely mechanisms for enhanced scattering of propa- 
gating vibrational excitations in a-RE^Sii-^ compared to 
a-Si: effects of the columnar structural features, and rat- 
tling of the caged Y or Gd dopants. The measurements of 
C and k presented here indicate that both mechanisms 
play a role in the thermal properties of the rare-earth 
doped amorphous silicon. 

Figure EJi compares the estimated mean free path of 
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FIG. 6: Estimated mean free path for vibrations vs. T. a) 
a-Yi2Si88 (heavy solid line) compared to a a-Si film (dashed 
line) grown by the same technique. Solid lines indicate es- 
timated T-dependence. Below ~ 15 K, where Xdom is ~ to 
the spacing of the columnar features in the doped materials, 
the temperature dependence of the MFP in a-Si increases, b) 
Filled (La .75Tho.2Fe 3 CoSbi2) and unfilled (CoSb 3 ) skutteru- 
dite antimonides. 27 Filling the skutterudite cage structures 
with La reduces the estimated MFP over the whole measured 
temperature range. 



vibrations for amorphous Si and a-Y^Sii-^ films. This 
is given by 3k/CoVD, where vp is the estimated Debye 
velocity and Co is the corresponding specific heat con- 
tribution from propagating modes, estimated here using 
D = 487 K measured for a-SiJA In a-Yi 2 Si 88 MFP is 
roughly oc T throughout the measured temperature 
range, while in a-Si the exponent increases below ^ 15 K 
so that MFP oc T~ 2 08 . 15 K is the temperature where 
^dom is roughly equal to the spacing of the columnar 
features seen in Figs. ^> and |2p. Fig. [BJd makes a simi- 
lar comparison between filled and unfilled skutterudites. 
The mean free path is again estimated using Sk/CoVD, 
where k, Vd, and Cd values reported by Sales et alS Fill- 
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FIG. 7: C'/Cd plot compares contributions to C in excess 
of the Debye model for crystalline Si^- a-Si^" and a-YgSigi. 
The low C/Cd for a-Si compared to crystalline Si is atypical 
for amorphous systems, where the peak in C/T 3 is typically 
larger than in the corresponding crystal. The larger C/Cd 
in a-YgSigi is closer to expected behavior for amorphous sys- 
tems. 



ing the cage in the skutterudite crystal with the heavy 
La dopant reduces the estimated mean free path over 
the whole measured temperature range. One interpreta- 
tion of our data is that the overall suppression of k in 
a-Y^Sii-a; is due to the effect of the rattling of caged 
rare-earth dopants, while the reduced temperature de- 
pendence of the MFP below 15 K is the effect of the 
columnar structural features. This structural scattering 
need not be directly analogous to grain-boundary scat- 
tering in a crystal, but rather the effect of a softer region 
of the matrix occurring due to density fluctuations with 
the 3 — 4 nm spacing seen in Figs.^ 3 and 03. 

Figure [7| compares the ratio of the measured C to the 
Debye contribution, Cd, for crystalline silicon, a-Si, and 
a-YgSigi. In typical amorphous solids such as a-Si02, C 
is significantly larger than Cn, as a result of the excess 
modesi We previously reported that the situation is very 
different in a-Si, where C more closely matches the Debye 
model than does crystalline Siii& As shown in Fig.0 the 
additional non-propagating vibrational states introduced 
by adding Y to the matrix cause C/Cd for a-YgSigi to 
slightly exceed the value for crystalline Si, approaching 
the expectation for a typical amorphous material. It is 
also interesting to note that for a-YgSigi, a-Yi3Sig7, and 
a-Y2iSi7g, the height of the peak in C/T 3 , P c , scales with 
its position in temperature, T max , as P c oc T~^, which 
agrees with the scaling observed by Liu and Lohneysen 
for a wide range of amorphous solids 4& Our recent a-Si 
measurement does not match this scaling behavior par- 
ticularly well, indicating again that a-Si is a somewhat 
atypical amorphous material, while the additional non- 
propagating vibrational states cause a-Y^Sii-a; to show 



more typical behavior. 

The exact nature of these states is unknown, though 
it seems likely that excess modes are added due to lo- 
cally softened Si cages around the heavy dopants or to 
anharmonic rattling of the dopant in the cage or both. 
It is clear from the reduced k that these excess modes do 
not carry heat. Several authors have reported a similar 
increase in the specific heat of filled skutterudites when 
compared to the "empty" host crystal. In the case of 
the skutterudites, the excess C can be explained by the 
contribution of one or more Einstein modes, which each 
contribute a term C E = (9 E /T) 2 e^/ T ^ / (e^/ 7 ) - l) 2 , 
where #e is the Einstein temperature. The observa- 
tion of these Einstein contributions provides direct ev- 
idence of the localized "rattling" of the filler atom, and 
is corroborated by inelastic neutron scattering data, res- 
onant ultrasound spectroscopy, and simulations. 30 32 - 33 
Following the "rattling" atom analogy, we fit the mea- 
sured specific heat of a-YgSigi to the equation C = C(a- 
Si) + J 4iCEi + J 4 2 CE2 + ^l3C'E3- A fit of this type assumes 
that the filling of the cages has little effect on the elastic 
properties of the host material. Though we believe this 
is true to good approximation for the a-YgSigi sample, 
available evidence suggests a possible reduction in #d for 
more heavily Y-doped a-Sii 3 ^ This type of effect can have 
important implications for the physical interpretation of 
the resulting fit parameters^ The solid line in Fig. 0] 
represents the fit with Ay = 1.0 J/mol K, 6>ei = 194 K, 
A 2 = 0.83 J/mol K, 6> E2 = 100 K, and A 3 = 0.086 J/mol 
K, 6*E3 = 52 K. The data can be fit somewhat less well by 
two Einstein modes (~ 108 and 51 K) , but is very poorly 
modeled with a single Einstein contribution^ As a com- 
parison, Keppens, et al. used Einstein modes at 70 and 
200 K to explain the contribution of the rattling La atom 
in Lao.gFe3CoSbi2r^ while Hermann, et al. needed only 
a single Einstein mode at 53 K to describe thallium rat- 
tling in similar antimony skutteruditesi 3 ^ The similarity 
of these results to the Y-doped amorphous Si data pre- 
sented here suggests that similar physics drives the excess 
modes and reduced k in these two classes of materials. 

Regardless of the exact nature of the contributions to 
C in a-Y^Sii-z, the combined effect of the addition of ex- 
cess modes and the enhanced scattering is an amorphous 
material with a large bump in C/T 3 but no correspond- 
ing plateau in k. We conclude that these two phenomena, 
which are often both observed in a given amorphous ma- 
terial, are not necessarily the work of the same physical 
mechanism. This suggests an interesting future experi- 
ment. It is currently a matter of debate whether the k 
plateau and C/T 3 bump and the related excess modes 
can be explained within a single theoretical framework 
with the TLS that dominate the low T properties of most 
amorphous insulators. Our work shows that addition of 
heavy Y dopants to a-Si, a material which has few ex- 
cess modes and very low contributions from TLS, adds 
excess modes to the material. Study of a-Y^Sii-a; below 
2 K could indicate whether TLS have returned with the 
addition of Y, which would provide evidence of correla- 
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tion between TLS and excess modes. This also suggests 
that addition of heavy dopants could provide a poten- 
tially tunable method for adding excess modes to a-Si, 
allowing systematic study of the vibrational excitations 
in amorphous insulators. 

IV. CONCLUSIONS 

In summary, we measured specific heat and thermal 
conductivity of Gd- and Y-doped amorphous Si thin 
films. Addition of the heavy dopant atoms introduces 
scattering which reduces the thermal conductivity below 
values measured for pure a-Si. At the same time, a large 
number of excess vibrational modes are added, result- 
ing in a bump in C /T 3 which is much larger than that 
in the pure amorphous material and also larger than in 
crystalline Si, in better agreement with expected behav- 
ior of C in amorphous insulators. This bump can similar 
manner as the contribution of rattling modes to filled 



antimony skutterudites, suggesting that similar physics 
drives the thermal behavior of these rather different sys- 
tems. Furthermore, the ability to add excess modes to 
a-Si suggests that continuing study of heavy-atom doped 
amorphous silicon could enable systematic probe of the 
correlation between tunneling systems and excess modes 
in amorphous materials. 
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